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Abstract 

We consider the many-particle quantum mechanics of anyons, i.e. 
identical particles in two space dimensions with a continuous statistics 
parameter a £ [0, 1] ranging from bosons (a = 0) to fermions (a = 1). 
We prove a (magnetic) Hardy inequality for anyons, which in the case 
that a is an odd numerator fraction implies a local exclusion principle 
for the kinetic energy of such anyons. From this result, and motivated 
by Dyson and Lenard's original approach to the stability of fermionic 
matter in three dimensions, we prove a Lieb-Thirring inequality for 
these types of anyons. 

1 Introduction 

The concept of identical particles and associated particle statistics lies at 
the foundations of quantum mechanics. It arises as a consequence of the 
non-observability of particle interchange and the fact that states in quan- 
tum mechanics are represented by rays in a complex Hilbert space, i.e. only 
determined up to a complex phase. A quantum mechanical state describing 
N distinguishable particles^ moving in M. d is represented by an N -particle 
wavefunction, i.e. a square-integrable complex- valued function u £ L 2 (M. dN ) 
defined on N copies of M. d , or equivalently, by an element of the tensor prod- 
uct space L 2 (M. d ) derived from the one-particle Hilbert space L 2 (M. d ). 
Upon restricting to identical particles, the freedom of choice of particle 
statistics stems from the fact that only the amplitude of the wave- 

function — describing (the square root of) the probability density for mea- 
suring the specific configuration x = (x\, X2, ■ ■ ■ ,xn) of particle positions 
xj £ M d — is observable, but not the exact phase u(x)/\u(x)\. Hence, since 
there should be no observable difference between the particle configuration 



"This work was partially supported by the Danish Council for Independent Research. 
1 We will for simplicity restrict to scalar non-relativistic particles, i.e. point particles 
without internal symmetries and spin. 
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x = (. . . , Xj, . . . , Xk, ■ ■ ■) and the one x' = (. . . , x^, ■ ■ ■ , xj, . . .), with parti- 
cles j and k interchanged, the amplitudes must be the same, but the phase 
may differ, as expressed by 

u ( x ') = e iaw u{x), (1) 

with a £ [0,2). In three and higher dimensions one finds that the only 
two possibilites for such a phase are e ia7T = ±1, corresponding to bosons 
(such as photons) with the plus sign, and fermions (such as electrons) with 
the minus sign. It is in fact sufficient to consider the permutation group 
Sn acting on the TV-particle Hilbert space L 2 (M. dN ), and the wavefunctions 
describing identical bosons respectively fermions are then given by the com- 
pletely symmetric resp. antisymmetric TV-particle wavefunctions, which in 
the latter case can be represented by the Hilbert space f\ N L 2 (M. d ). 

However, in two dimensions the statistics parameter a can be taken to 
be any real number in the interval [0,2) (or (—1, 1], by periodicity). Again, 
bosons correspond to a = and fermions to a = 1, while for a general choice 
of a the corresponding particles are simply called anyons (or, historically, 
just "particles obeying intermediate or fractional statistics" [3 El [2D]). This 
discrepancy between two and higher spatial dimensions is directly related to 
the fact that a punctured plane is not simply connected, while M. d \{0} is, 
for d > 3. Hence, for particles confined to the plane one has to consider con- 
tinuous interchanges of particles, forming a loop in the configuration space 
(together with the possibility of enclosing other particles in that loop), and 
the permutation group symmetry is replaced by the braid group Bjy, whose 
one-dimensional unitary representations determine the choice of statistics 
for identical anyons through the phase e ian . For a convenient and rigorous 
treatment of anyons, one can model them as identical bosons or fermions 
in the plane, but with a magnetic interaction of Aharonov-Bohm type be- 
tween each pair of particles, giving rise to the correct statistics phase as 
the particles encircle eachother. The quantum mechanical momenta of the 
particles, which for bosons and fermions are simply given by the gradients 
Vjii w.r.t. the particle positions xj, will then be replaced by covariant 
(magnetic) derivatives Dju (see (fT2|) below). 

Fermions in any dimension are special, since they satisfy the so-called 
Pauli exclusion principle. Namely, because of the antisymmetry of the wave- 
function, no two particles can occupy the exact same state, expressed simply 
with the help of the wedge product as uq A uq = for any one-particle state 
uq G L 2 (R d ). An important and non-trivial consequence of this is the cele- 
brated Lieb-Thirring inequality, which given a scalar potential V on ~M. d can 
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be summarized as 



V / (\Vju\ 2 + V(x j )\u\ 2 ) dx 
~i JR dN 

N-l 

> -El^WI ^ -Cf d \V-(x)\ l +i dx, (2) 

where A&(/j) denote the negative eigenvalues (ordered with decreasing mag- 
nitude) of the one-particle Schrodinger operator h := — A R d + V(x) acting in 
L 2 (R d ). The first inequality just expresses the fact that, because of the Pauli 
principle for fermions, the lowest possible energy (l.h.s. of ([2])) is obtained 
when the particles assume the states corresponding to the lowest N eigen- 
values Afc(/i), i.e. when u = f\k=o Uk is an antisymmetrized product of those 
one-particle eigenfunctions Uk{h). The second inequality holds uniformly in 
N and concerns the trace over the negative spectrum of the one-particle op- 
erator h. Ever since the first proof of ([2]) in 1975 by Lieb and Thirring |12j . 
who used this result for a simplified proof of stability of fermionic matter 
(see also there has been a lot of activity in the mathematical com- 

munity aiming to generalize this type of spectral estimate for one-particle 
operators in various directions. Furthermore, the Lieb-Thirring inequality 
([2]) (disregarding the intermediate sum over eigenvalues) is equivalent to the 
kinetic energy inequality 

N 

T:=V / \V jU \ 2 dx > C$ / P {x) l+1 ddx, (3) 
~[ Jm dN Jm. d 

which can be interpreted as a strong form of the uncertainty principle for 
fermions, because of the way it bounds the total kinetic energy T in terms 
of the one-particle density 

N 

p(x):=y^ \u(xi,...,Xj = x,..., x N)\ 2 T\dx k (4) 

3=1 k^j 

of the wavefunction (always assumed to be normalized to f Rd N \u\ 2 dx = 1). 

Bosons, on the other hand, do not satisfy any exclusion principle. They 
can all be put in the same state, e.g. u = uq <S> ■ ■ ■ <S> uo, and hence cannot 
be expected to satisfy the inequalities ([2]) or Q. In fact, the best we can 
do is to treat them as N copies of a single particle satisfying ([2]), and hence 
the above inequalities hold only in the weaker form 

at 

^2 / {Nju\ 2 + V( Xj )\u\ 2 ) dx > -C% T N I \V„{x)\ l+ i dx, 
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resp. (cp. Appendix B) 



" r nK , 



E 



which now only encodes the uncertainty principle, without the extra gain in 
([3]) due to statistics. These bosonic inequalities become trivial as N — > oo. 

In contrast, not much has been known about the spectral and statistical 
properties of many anyons for < \a\ < 1. Not even the ground state energy 
of an otherwise non-interacting gas of anyons has been rigorously estimated^. 
The main difficulty lies in the fact that many-anyon wavefunctions cannot be 
simply related to one-particle wavefunctions in the same way as for bosons 
and fermions. As will be seen explicitly below, the Hamiltonian operator 
Hq = X^'=i Dj describing the kinetic energy of N free anyons is not just a 
free Laplacian acting on totally symmetric or antisymmetric wavefunctions, 
but involves long-range magnetic interactions between all the particles. In 
particular, we cannot reduce our study to the relatively simple case of a one- 
particle Schrodinger operator h. The aim of the present paper is to address 
this gap in knowledge concerning intermediate anyon statistics, as well as 
the current lack of techniques to study the spectral theory of such quantum 
mechanical systems. 

Our first main result concerns a magnetic many-particle Hardy inequality 
for N anyons, which when considered on the full two-dimensional plane M 2 
reads: 

V I \D jU \ 2 dx > V f , H ' {2 dx, (5) 

f^jR 2N N Jwl 2N \ x i ~ x j\ 

with the statistics-dependent constant 

C a , N ■= min min|(2p + l)a-2g|. (6) 

p=o,i,...,iV— 2 gez 

A stronger form of ([5]), given as Theorem 2] below, and valid for any convex 
subdomain C M 2 , is shown to produce a local form of Pauli's exclusion 
principle for anyons (given as Lemma [7] below), whose strength depends on 
the large- iV behavior of the constant C a> jv- Although C&,N is clearly non- 
zero for all N whenever a is irrational, we find that infTveN Ca,^ = 0, unless 
a = with /i and v relatively prime integers and [i odd, in which case we 
shall see that infjv g ^ C a ^ = \. 

For such odd numerator fractions a = ^ , the energy given by the local 
Pauli exclusion principle for arbitrary numbers of particles is of a similar 
form as for fermions, but with an extra factor \ (depending wildly on the 



2 Note that this is trivial in the case of bosons, and a simple exercise in the case of 
fermions. The problem for anyons has been attacked by many authors through various 
approximations (see the references in the books and reviews [H 171 HOl 1161 121] ). 
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statistics parameter). Our second main result is that this local bound is 
sufficient to produce a Lieb-Thirring inequality for this type of anyons. In 
our approach we have been inspired by Dyson and Lenard's original proof of 
the stability of fermionic matter in three dimensions [2] (from 1967, before 
the advent of the Lieb-Thirring inequality), in which the only place where 
the Pauli principle came in was through such a local bound for the energy. 

Theorem 1. For a normalized wavefunction u of N anyons on M 2 , with 
odd-fractional statistics parameter a = ^ (i.e. a reduced fraction with p 
odd), we have the kinetic energy inequality 

N 

T:=V / \D jU \ 2 dx > C K ^ [ p(x) 2 dx, (7) 

fr[ Jr 2n » Jm. 2 

where p is the corresponding one-particle density and hence the Lieb- 
Thirring inequality 



N 

E / 



\D jU \ 2 + V(xj)\u\ 2 ) dx > -C LT v 2 [ \V-(x)\ 2 dx, (8) 

Jr 2 



for any real-valued potential V on M 2 . Here Ck and Clt are universal 
positive constants that can be given explicitly. 

In particular, this implies that the total kinetic energy T per unit area 
for a non-interacting (apart from the statistical interaction) gas of anyons 
with odd- fractional statistics a = confined to an area A, is bounded 
below by 

T p 2 

-7 > Ck-o, 

A u z 

where p := N/A is the average density of the gas. Other implications of this 
Lieb-Thirring inequality, such as for interacting anyon gases, as well as the 
sharpness of this result concerning even-fractional and irrational statistics 
and its physical interpretation, will be discussed in a forthcoming paper |14j . 

The structure of this paper is as follows. In Section 2 we fix the notation 
and briefly recall the general theory of particle statistics. Our fundamental 
many-particle Hardy inequality for anyons is proven in Section 3 based on 
a pairwise relative parameterization of the configuration space, combined 
with a local magnetic Hardy inequality which takes into account the under- 
lying symmetry between the particles. Section 4 concerns the local gain in 
energy following from these Hardy inequalities, which we refer to as a local 
Pauli exclusion principle for anyons due to its direct similarity with the cor- 
responding local gain in energy for fermions. In Section 5 we use this local 
gain to prove the Lieb-Thirring inequality for anyons. In the appendices we 
have placed some suggestions for improvements of the local energy, as well 
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as a proof of a Lieb-Thirring inequality on cubes with Neumann boundary 
conditions. 

We emphasize that, although a Lieb-Thirring- type inequality for anyons 
could perhaps have been anticipated based on physical grounds (at least 
for some values of the statistics parameter a), this is from a purely mathe- 
matical perspective a highly non-trivial extension of the usual Lieb-Thirring 
inequality ([2]) since the relevant operator 

N 

H = Y J {D 3 (xf + V{x 3 )), 

3=1 

with Dj{x) being the differential operators given explicitly in (|12p below, is 
now a strongly interacting magnetic many-particle Hamiltonian. 

Acknowledgements. D.L. would like to thank Oscar Andersson Forsman for 
discussions related to Appendix A. 

2 Preliminaries 

2.1 Particle statistics in two dimensions 

In this subsection, which is not a prerequisite for understanding the rest 
of the paper, we give a very brief recap of the general theory of particle 
statistics in d > 2 dimensions. For more details we refer to the original 
reference [9], the review articles [HQS], and the books QUI [21] on anyons. 

The classical configuration space of N identical particles in M. d is formally 
given by 

X? :={R dN \B) /S N , 

where we have excluded all coincidences of the particles, i.e. the diagonals 

B := {x G R dN : Xj = x k for some j ^ k}, 

and the symmetric group Sn acts on the N copies of M d in the obvi- 
ous way. (The center-of-mass coordinate X := Ylj x j can be trivially 
factored out, = M. d xl* rel , leaving the relative configuration space 
X^ vel = ({x G R dN : X = 0} \ B) /S N .) For d > 3 the fundamental group 
of X^ 1 is Tti{X^) = Sn, whereas for d = 2 it is the braid group on N 
strands, tti^X^) = B^. Wavefunctions of ./V identical particles are defined 
as square-integrable complex-valued functions on X^ with appropriate glu- 
ing conditions (recall the physical requirement ([I])). Hence, these can be 
viewed as sections of a complex line bundle over X^ . There are natural flat 
connections on such line bundles, taking the trivial connection locally on 
M. dN (note that the parallel transports could be globally non-trivial as there 
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e i(2p+l)a-K 



Figure 1: One- resp. two-particle interchange loops with corresponding braid 
diagrams (where we can think of time as running upwards) and phases. 

are non-trivial transition functions between local regions), and every such 
connection defines a unitary one-dimensional representation of the funda- 
mental group iri(X^). For d > 3 there are only two such representations, 
the trivial one corresponding to bosons, and the sign on Sn corresponding 
to fermions. For d = 2 the unitary one-dimensional representations of the 
braid group are parameterized by a real number a G [0,2), where every 
generator in Bjy corresponding to a counter-clockwise interchange of two 
neighbouring strands is represented by the phase e ia7T . In particular, this 
will imply that 



corresponding to continuous counter-clockwise interchange of two particles 
Xj and Xk, with the interchange loop enclosing precisely < p < N — 2 
other particles x.^, . . . ,Xi p . On the other hand, if a single particle Xj is 
taken along a simple loop which encloses p other particles, then a phase 
factor e l2 'P a7T will be picked up. See Figure Q] for these examples, and for a 
glimpse of how the exact phases arise through braid diagrams, with each 
elementary counter-clockwise braid contributing a phase e ia7r . If a = or 
a = 1 the representations depend only on the permutations of the particles 
and not on the braid, and we are back to the case of bosons and fermions. 

We shall denote the Hilbert space of iV-particle wavefunctions with a 
given statistics parameter a by . In the case of bosons or fermions we 
may identify such sections of line bundles with complex-valued functions on 
M. dN that are either totally symmetric or totally antisymmetric, and hence 
we obtain the usual spaces of iV-boson, resp. iV-fermion wavefunctions. 

In general, the line bundles corresponding to different a are topologi- 
cally equivalent, but not geometrically equivalent. The map (singular gauge 
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transformation) 



u(x) ^ Yle^ a - a ^u(x), ^ : =arg^ L (l ) 

(where we pick an arbitrary real axis to identify R 2 with C) maps the line 
bundle corresponding to Qo to the line bundle corresponding to a. The 
natural flat connection on the ao-bundle will then not be mapped to the 
natural flat connection on the a-bundle, but gives rise to the non-trivial 
gauge potential 

Aj := -i /^e^-" ^ V Xj (j^e^-^A . (11) 

\k<i J \k<i / 

We can then think of the ao-bundle as providing a reference statistic, and 
we will choose ao to be either or 1 in order to model general statistics a 
in terms of bosonic or fermionic wavefunctions. In this way we can model 
iV-anyon wavefunctions as either totally symmetric or antisymmetric N- 
particle wavefunctions u G H^o — L 2 (M? N ) with covariant derivatives 

D 3 = iV, ■ A,. 

The advantage of taking this viewpoint is that we then can work exclu- 
sively on the redundant but uncomplicated configuration space IR 2 ^. The 
disadvantage, of course, is that we have to deal with a long-range magnetic 
interaction potential Aj. In the physics literature the former viewpoint is 
often referred to as the anyon gauge, while the latter is called the magnetic 
gauge. Although we will mostly stick to the magnetic gauge (with ao := 0), 
it is for the purpose of intuition very useful to have both pictures in mind. 

2.2 Notation 

In order to make precise our notation, we will denote by @a Q a the space 
of finite kinetic energy wavefunctions of ./V anyons with reference statistic 
ao £ {0, 1} and actual statistics parameter a, i.e. the totally symmet- 
ric/antisymmetric functions u E Ha — L 2 (R 2N ), where each particle is 
interacting with all the other through Aharonov-Bohm magnetic potentials 
of strength a — ao, in the following precise sense. The covariant derivative 
acting w.r.t. particle Xj (following from (|10p - (|lip ) is 

Dj := -iVj + Aj(xj) := -iV Xj + (a - a ) ^(xj - x^I, (12) 

where a 1 := a/\a\ 2 and a i— )■ al denotes counter-clockwise rotation of the 
vector a £ I 2 by an angle ir/2. We consider the semi-bounded quadratic 
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form 



u h-> q(u) := 

3=1 



Djit\ dx = 




(13) 



defined initially on u G C$°(R 2N \D) n %^ , i.e. the smooth totally sym- 
metric/antisymmetric and square-integrable functions on M 2N with support 
away from diagonals (this amounts to what is in the literature sometimes 
referred to as a hard-core condition on anyons). The space of finite kinetic 
energy iV-anyon wavefunctions Ql^ a is then defined as the domain of the 
closure of this quadratic form q on rl^ . Note that there is also an associated 
smaller space a defined as the domain of the Priedrichs extension, i.e. 
of the self-adjoint operator Hq = ^ - D 2 on T~L^ C L 2 (M. 2N ) associated to 
the closure of the quadratic form (|13|) . and that for ckq = a = 0, 



where H k denote the Sobolev spaces of k partial derivatives in I? . In our 
proofs we will always use the denseness of the smooth functions with compact 
support in these spaces, and pick such representatives without taking explicit 
limits. 

In the following, open resp. closed balls of radius r at a point x will 
be denoted B r (x), resp. B r (x), and the characteristic function of a set A 
is denoted xa- Given a real- valued function or expression /, we define the 
non-negative quantities f± := max{0, ±/}. 

3 Hardy inequalities for anyons 

We can mention [6], where a many-particle Hardy inequality has been de- 
rived for anyons (see Theorem 2.7 in [6]), but which is unfortunately not 
sufficient for our purposes^. The weakness stems from the fact that, so far, 
only single-particle movements have been taken into account, which only 
captures some of the symmetries involved (cp. Figure [T]). In order to ar- 
rive at something non-trivial, at a minimum for the special case a = 1 of 
fermions, we need to consider a relative Hardy inequality (cp. e.g. Lemma 
4.6 in [6] which one could view as a relative Hardy inequality for a pair 

3 Note, for instance, that the corresponding Hardy constant in [6], 



is zero for a = 1 (and any a € Q for TV large enough), and that in any case Dn.c* < iV -2 . 
Hardy inequalities for interactions of anyonic type have also been considered in [15] . 
but these are of single-particle type and hence also do not take the underlying symmetry 
between particles into account, as well as have an unclear dependence of the corresponding 
constants on the positions of the Aharonov-Bohm fluxes. 




H 2 (m 2N )nn^ c H 1 (m. 2N )nn£ = ®t 
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of fermions in any dimension). Also crucial for our approach is the fol- 
lowing extension of a class of well-known two-dimensional magnetic Hardy 
inequalities (see [HI HJ [15] ) , where we take the underlying symmetry of the 
wavefunction and gauge potential into account. 

Lemma 2 (Magnetic Hardy inequality with symmetry). Let 17 = Br 2 (0) \ 
Br 1 (0), R2 > Ri > 0, be an annular domain in M. 2 , and let there be a 
magnetic flux inside -6^(0), determined on £1 by a vector potential a : 
17 — > M? , s.t. V A a = on 17 and f r a ■ dr = & for any simple loop T in Q 
enclosing i?/? 1 (0). Furthermore, assume that a is antipodal-antisymmetric, 
i.e. a(—r) = — a(r) for all r G 17, and let v S C°°(17) be a function on 17 
with antipodal symmetry, v(—r) = v(r) for all r 6 17. Then 

2 



\D r v\ dr > min 



2k 

2ir 



7-79 dr, 



where D r := — iV r + a(r). 

Alternatively, if v is antipodal-antisymmetric, v 
r 6 17, then 

2 



|D r i;| dr > min 



dr. 



(14) 



-v(r) for all 



(15) 



Proof. There exists a gauge transformation v ^ v = e %x v such that |-D r ^| 2 = 
|(— iV r + a)v\ 2 , where a 1— > a(r) := §^r~ l I on 17. Note that x{ r )i being 
the integral of a difference of two gauge potentials a(r) and a(r), both 
antisymmetric w.r.t. r 1— > — r, must be symmetric under this antipodal 
map. Hence, if v is antipodal- (anti)symmetric, then so is v. 

Now, we write the gauge-transformed l.h.s. of (|14p in terms of polar 
coordinates (r, (p), 



[ \D r v\ 2 dr = [ 
Jn Jo 



2tt r R 2 



\d r v\ 2 r drdip- 



2tt r R 2 



Hi 



1 



-idu + 



-) 

2ir J 



r drdip. 



Considering only the last term involving cL>, and Fourier expanding v(r) on 
17, 

v(r, ip) = v{r = rexe* 1 ) = -L ^ v k {r)e i2kip 



(note that we have here used the fact that v(r,ip + n) = v(r, (p) for all ip), 
we find 

2 r I |2 



R 2 
Ki 



2A; + — 
2vr 



2 l~ 12 



r <ir > min 

fcez 



2fc 

2vr 



dr, 



n 



and hence arrive at the inequality (|14|) . For the case of antipodal-anti- 
symmetric v, we can Fourier expand v in odd powers of e ikip and arrive at 
(USD. □ 
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Now, let us for simplicity first apply this to the case of only two anyons 
and prove a relative Hardy inequality for this system. 

Lemma 3 (Relative two-anyon Hardy). Let fi be an open convex set in M 2 
and let u G S>Q a be a two-anyon wavefunction. Then 



where 



\Dou\ ) dx\dx 2 > 2 min la 

' k£Z 



2k\' 



u 



non Fi - «2 1 



dxidx 2 , 
(16) 



fiofi := {(£Ci,£C 2 ) G fi 2 : ||sci - x 2 \ < dist(±(a5i + x 2 ), fi c )} (17) 
(m particular, fi o fi C fi 2 , unless fi = M 2 J. 

Proof. Let us introduce the center-of-mass -R := i(aci + 2:2) and the relative 
coordinate r := ^(x\ — x 2 ). Furthermore, let v(R;r) := u(R+r, R—r), and 
observe that the bosonic symmetry of u G T-Lq implies v(R;—r) = v(R;r) 
for all R G fi and r G M 2 s.t. < \r\ < dist( J R,fi c ) =: S(R) (possibly 
infinite). Then, by V R = Vi + V 2 , V r = Vi - V 2 , the l.h.s. of CED equals 




nJB S{R) (o) 



a 



+ 



Q 



4 /U 



2drdR 



J 



\V R v\ 2 drdR- 



-iV r + ar~ v\ drdR. 



For the last integral w.r.t. r we can then apply Lemma [2] (where no gauge 
transformation is necessary in this case), and thus find 

I {\Diu\ 2 + \D 2 u\ 2 ) dxidx 2 > min I a - 2k\ 2 [ [ •pL drdR, 



which implies (|16p . 



□ 



An equivalent way to think about this result in terms of the anyon gauge, 
i.e. the space % 2 , is that we can make a symmetric interchange of the two 
anyons by rotating them around their common center-of-mass. After half a 
turn — already then completing a loop in the relative configuration space 
^2 vei = \{0})/ r «._ r — we can compare the function values, and the 
condition ([9]) on parallel transport tells us that we should pick up a phase 
e tan . We are therefore given a function v(r) on, say, the upper half-plane 



with a periodic boundary condition v[—rei) 
Hardy inequality follows. 



e ta7T v(rei), and hence the 
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This is straightforwardly extended using Lemma [2] to cases when there 
are also fluxes inside the interchange loop. Intuitively, as we encircle p 
anyons with a symmetric two-anyon interchange loop we pick up a phase 
e i(2p+i)air anc j h ence a corresponding Hardy inequality. More precisely, this 
results in the following relative many-particle Hardy inequality for anyons. 

Theorem 4 (Many-anyon Hardy). Let Q be an open convex set in M 2 and 
let u S ^g a be an N-anyon wavefunction. Then 

I ir\Dju\ 2 dx> — [ y\ — — — n Xnon(xi,Xj)dx, (18) 
Jn*^, N Jqn^. \xi-Xj\ 2 

3=1 *<J 
with Q o $7 defined in (|17p . and 

C Q , N := min mm |(2p + l)a - 2g|. (19) 

p=0,l,...,N— 2 q£Z 
Proof. We use that, for any z = (zj) G C dN , 



i=i l<j<k<N 



N 
3=1 



(20) 



Applying this identity with Zj := Dju(x) and dropping the second term, we 
see that the l.h.s. of ()18j) is bounded below by a sum of (^) integrals of the 
form 

/ / \(Dj ~ D k )u\ 2 dxjdx k \\ dx h (21) 

Jn N - 2 J{xj,x k )&Q? 1=1 N 

Mi,* 

Hence, we consider for each fixed choice of the N — 2 variables (xi) the 
remaining configurations of the pair (xj, x^) on Vl o C which we again 
parameterize by a center-of-mass R := |(a;j + aJfc) and relative coordinate 
r := \{xj — x k ). For each _R 6 SI we can then split up the final param- 
eterization of r G Bgr R \(0) \ {0} into annuli, the first annulus extending 
from r = 5q := to r = Si, defined to be the distance from R (where there 
could possibly be one particle x{) to the next closest particle xy (or possibly 
several particles situated on the same distance from R). The next annulus 
then extends from r = 8\ to the next greater distance 82 from R to any par- 
ticles, and so on, until we reach the boundary of the domain (or infinity) at 
r = 8m '■= 8(R). Now, on each annulus A m := B$ m \ Bg m _ x , m = 1, . . . , M, 
we have 

(Dj - D k )u = (-iV r + ar~ l I + a(R; r)) v, 

where v(R;r) := u(R + r,R — r) is antipodal-symmetric in r, while the 
gauge potential a(R; r) := A{R+r) — A{R—r) is antipodal-antisymmetric, 
with A[x) := a YLi^j k( x ~ x i)^ 1 ^ being the magnetic potential at x £ 
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Ill > 



R + A m from all the other N — 2 particles X[ . Note that V r A a = on A 
and that f r a- dr = 4irap m for any simple loop T C ^4 m enclosing Bg m _ 1 , 
where p m is the number of particles x\ inside We can therefore apply 

Lemma[2l with the total flux $ in the disk Bg ml given by $ = 2ira(l+2p m ). 
Hence, 

/ \{Dj - D k )u\ 2 dr > 4mm\(2 Pm + l)a - 2q\ 2 [ ■ ^ ^ dr 

JA m <?ez JA m \Xj-Xk\ 2 



and proceeding similarly for all annuli A m , all points i?, and all pairs (j, k), 
we obtain the inequality (fT8|) . □ 



Note that this theorem coincides with Lemma [3] for N = 2, and with 
Theorem 2.8 in [6] for fermions on M 2 (since C a= \^ = 1 for all N > 2). 
We also note that C a= o,Ar = is the optimal constant for bosons since 
(|18p concerns the Neumann form, i.e. we could in this case (and if Q has 
finite measure) take u to be constant so that the l.h.s. is identically zero. 
If one considers the Dirichlet form on the other hand, a non-trivial many- 
particle Hardy-type inequality for bosons or distinguishable particles in two 
dimensions (also with a constant ~ N^ 1 , although with logarithmic factors 
in the potentials) was derived in |13j . 

Concerning intermediate statistics, we have the following observation: 

Proposition 5. The infimum 

C a := inf C aN = inf \(2p + l)a - 2q\ 

NeN ' p,q& 

equals 



a 



a 



if a = £ with fi G Z, v E N+ relatively prime and fi odd, 
otherwise. 



Hence, C a is strictly positive whenever a is an odd numerator fraction, but 
zero otherwise. 

For an upper bound in the non-trivial cases we will use the following 
fact: 

Lemma 6. Given a, b £ Z coprime, where a is odd, there exist integers x, y 
such that x is odd, y is even, and ax + by = 1 . 

Proof. Since a, b are coprime, we can by Euclid's algorithm find integers x, y 
s.t. ax + by = 1. We then use that ax + by = a(x + kb) + b(y — ka) for 
any k S Z and consider the different possibilities. If b is odd then, since a is 
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0.8 



0.6 



0.4- 



0.2- 



♦* ♦ * ♦ * ♦ 



0.2 



0.4 0.6 
a 



0.8 



Figure 2: A sketch of the behavior of C a as a function of a. 



also odd, either x is odd and y even and we are done, or x is even and y is 
odd, but then we can choose k odd and so we are done. In the case that b 
is even, then either x is odd and y even and we are done, or both x and y 
are odd. In the latter case we can again choose k odd and we are done. □ 



Proof of Proposition^ Assume that a = with fi € Z, v 6 
and consider first the case of odd numerators: \i = 2k + 1, fc£Z. 



coprime, 
Then 



1 



1 



|(2p + l)a - 2q\ = -\(2p+ l)(2fc + 1) - 2gi/| > -, Vp,q£Z 



i.e. C a > -, since the last absolute value is a difference between an odd and 
an even integer. Furthermore, Lemma [6] guarantees the existence of p, q G Z 
such that this absolute value is equal to one, hence C a = ~. 

In the case of even numerator fractions, i.e. /i = 2k, k £ Z, we have 
\(2p + l)a — 2q\ = \\(2p + V)k — vq\. As v necessarily is odd, we can choose 
2p + 1 = v and q = k, and hence C a = 0. 

In the irrational case a E M\Q, we can find an infinite sequence of 
p,q £ 7L such that (see [19]) 



a 



2q 



2p+l 



< 



1 



(2p + 1) 



and hence C a = 0. 



□ 



In Figure [5] we have plotted a sketch of the dependence on a of the 
large- N constant C a . It should be emphasized that for small numbers of 
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particles (relative to a fixed a), the Hardy constant C a> N is typically non- 
zero. Its graph can be obtained by cutting out wedges of slopes 2p + 1, 
p = 0,1, . . . , N — 1, at every rational point with denominator 2p + 1 and 
an even numerator. In particular, C a ^N > for all irrational a, and for 
even numerator fractions a = ^ E (0, 1) the constant is strictly positive if 
and only if N < Let us also point out that we could just as well have 
chosen to work with the fermionic reference statistic, i.e. «o = 1> f° r which 
the corresponding statistics-dependent constant is 

C PtN := mm mm \(2p + - (2q + 1)|. (22) 

p=0,l,...,iV— 2 gg^ 

We have here denoted the strength of the statistical interaction by (3 := a—1, 
and the values of (3 for which this constant is bounded away from zero for 
all N are the fractions where either the numerator or denominator is even. 

4 A local Pauli exclusion principle for the kinetic 
energy 

Note that for odd- fractional statistics, i.e. for non-zero values of C a , the 
total constant in the many-anyon Hardy inequality (|18j) still tends to zero 
like N~K NaivehJ!, this would be insufficient for a non-trivial bound on the 
energy of the anyon gas in the thermodynamic limit, since the energy per 
area (L 2 , say, to which the gas is confined) due to Theorem 2] yields 

T 1 4C 2 (") f . l2 , N - 1 2 

T2 > I 2-^-^-yH 2 dx = const.^.p 2 ^0, 

as N — > oo, with fixed density p := N/L 2 . However, Theorem [3] is stronger 
than that, and we can choose not to consider this global gain in energy due 
to statistics directly, but rather its local implications upon cutting the space 
up and employing Neumann boundary conditions. This local approach is 
in the spirit of Dyson and Lenard's original proof of the stability of matter 
[2]. By choosing the size of such local regions appropriately, the energy can 
then be lifted to a stronger bound on the full domain R 2 (see (|3ip and (|36p 
below). 

Central for this approach is the following version of Lemma 5 in [2] for 
anyons. We refer to this as a local exclusion principle for the kinetic energy 
of anyons, since it implies that n > 2 anyons must have positive energy and 
therefore cannot all occupy the lowest zero energy state. 



4 Note that the sharp large-N behavior of this constant is not clear even in the fermionic 
case, cp. [6l|3]. We also remark that for d > 3 the corresponding sharp behavior is known 
[3], and does imply a non-trivial bound for the ground state energy of the fermion gas. 
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Lemma 7 (Local energy / Pauli principle). Let u G &Q a be a wavefunction 
of n anyons and let f2 C M 2 be either a disk or a square, with area Then 



[ J2\D jU \ 2 dx > (n-l)^p f \u\ 2 dx 



(23) 



where cn is a constant which satisfies cn > 0.477 for the disk, and cq > 0.358 
for the square. 

Proof. Note that by rescaling, we can in the following assume that Vt = B\ (0) 
or $7 = [— 1, l] 2 . We shall first consider the case of the disk, and then point 
out what needs to be changed for the square. 

Due to Q o O C Q 2 , the bound given by the many-anyon Hardy inequality 
(|18|) is unfortunately not sufficient as it stands, and we need to modify the 
approach in the proof of Theorem 0] to take the whole two-particle domain 
S7 2 into account. Instead of (|20p we shall therefore use 

n 



E n 2 = E ( ^zi (N 2 + i^i 2 ) + - Zk \ 2 ) + 1\ E^f' ( 24 ) 

i=i j<k v 7 i 



with < k < 1. The last term is again thrown away while the middle one 
is employed as in the proof of Theorem 0] to produce a Hardy potential in 
terms of relative coordinates r and R. However, for the first two terms we 
instead use that |-Djtt| 2 > Vj|u| (diamagnetic inequality). Ignoring k for 
a second, we are hence interested in the infimum of the ratio 

J 2 ^xM 2 + \^ X2 U \ 2 + -^XB S(R) (o)( r )\ u \ 2 ^ dxxdx 2 J J Ju\ 2 dxxdx 2 

(25) 

over u £ H 1 ^ 2 ), which is certainly greater than the lowest eigenvalue of 
the Schrodinger operator (c := C Qjn 7^ in the following) 

H:=-A% 2 + f, f( Xl ,x 2 ):=c 2 g(\R\,\r\), (26) 

on Q 2 defined with Neumann boundary conditions on dQ, 2 , with 



g(R,r) :-- 



5- 2 (l-R)- 2 , R < R, r < 6(1 - R), 
r" 2 , R<R, 6(1 - R) < r < 1 - R, 

k 0, R > R, 



for some fixed cut-off parameters < 6, R < 1, to be optimized over later. 

Now, denoting by P the projection onto the constant function uq(x) := 
|J7 2 |~'I = 7T _1 , and Q := 1 — P, we have (— A^)P = 0, and for the first 
non-zero Neumann eigenvalue Ai = Xi(-A^), 

(-A&)Q > Ax(-A^) Q = A x (-A# ) Q = £ 2 Q, 



16 



where £ ~ 1.8412 denotes the first zero of the derivative of the Bessel func- 
tion J\. Furthermore, we have since 

(u, (PfQ + QfP)u) = (pPu, pQu) + (f^Qu, pPu) 
<2\\f?Pu\\\\fsQu\\ <fi\\f 1 2Puf+-\\f l 2Qu\\ 2 = (u^fiPfP+^QfQ)^, 

for u G L 2 (0 2 ) and fi > 0, that 

/ = (P + Q)f(P + Q) > (1 - ^)P/P + (1 - [i~ 1 )QfQ. 

2 

These operators are estimated according to ||Q/Q|| < ||/||oo = ^7f ^2 1 an d 
\\PfP\\= [ fu dx 1 dx 2 = - [ c 2 g(R,r)2drdR 

JO? K 

T Jo \J0 JS(1-R) r J 

= W Q+ln^P 2 , 
and hence 

7J > (-A#JP + (— Aq_)Q + (1 - a«)P/P + (1 - /U-^Q/Q 



> (l-^)2vrc 2 (l + 21nr 1 )P 2 P+ I £ 2 - 0" 1 - 1) ) Q. 



5 2 (1-R)\ 

With the parameter k from (|24p reintroduced into (|25p and (|26p . we bound 



H > 2vrc 2 (l - M )(l + 21u5- 1 )P 2 kP + ^ 2 (1 - «) - Q 
> 0.304c 2 , 

where the lower bound was found by numerical optimization, with /x = 0.851, 
5 = 0.54899, R = 0.54396, k = 0.499, £ 2 > 3.389. Summing up, we have 

>— f [ ((1 - K,)(\Dju\ 2 + \D k u\ 2 ) + Kf(xj,x k )) dxjdx k dx' 

]<k 

- ( n ~ 1 )"T?T' / M***, ( 27 ) 

with c n > 0.304 • vr/2 > 0.477. 
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In the case that O = [— 1, l] 2 D -Bi(O), we can use the same / and g as 
above (extended by zero outside -Bi(O) 2 ), so that the bound on ||Q/Q|| is 
unchanged, ||P/P[| is multiplied by the area ratio f , and (-A# 2 )Q > ^-Q. 
We find 

H > —(1 - p)(l + 2\n5- l )R 2 nP+ ( ^(1 - «) - ^ ~ 1^ Q 

- 2 v wv ; ^ 4 v ; 5 2 (i-i?) 2 y 

> 0.179c 2 , 

with (j. = 0.8879, 5 = 0.5451, R = 0.531, « = 0.52, and hence (EH) holds 
with the constant c fi > 0.179 • 4/2 = 0.358. □ 

Remark. More general convex domains can be treated in a similar way using 
[17j . Also note that the bound (|23p for the disk holds with cq = 7r£ 2 ~ 10.65 
in the case that a = 1 (compare with Lemma 5 in [2]). For possible ways of 
improving the constant cq in the general case, see Appendix A. 

5 A Lieb-Thirring inequality for anyons 

Given a normalized A-anyon wavefunction u 6 71^ we will in the following 
denote by 

N 

p(x):=) / \u(xi, ... ,Xj-i, x, Xj+i, ... ,x N )\ 2 TT dx k 

^— ' / TO 2(iV-l) - LJ - 

i=i l/R k?y 

the one-particle density, such that J r2 p(x) dx = N . We start by reformu- 
lating Lemma [7] in terms of p. 

Lemma 8 (Local Pauli principle). Let u G be an N-anyon wavefunc- 
tion on M? and C M 2 a simply connected domain on which ()23p holds for 
some constant cq. Then 

n . c C 2 / r \ 

*> = E /„„ W*»<«*> * £ nrf (/»' " V ' (28) 

Proof. Using that 

N 

l = n{xn{x k ) + {l- X n{x k )))= Yl Hxn{x k )]J(l - Xn(x k )), 
fc=i ac{i,...,jv} fce^ fc^A 

the l.h.s. of ([2HD is 
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We now apply Lemma [7] to each term in the first summation above, which 
involves a partition A of the N particles into n := \A\ of them being inside 
the domain fi, while the remaining N — n residing outside, and therefore 
whose contributions to the magnetic potentials Aj £ a can be gauged away. 
Thus, we find 



> 



AC{1,...,JV} 

a 

a,N 



{\a\ - 1) + j j \ u \ 2 n dx k n dx k 

I \j u \ 2 E {\A\-l)\{xn{x k )\{(l-xn(x k ))dx. 



AC{1,...,N} keA k£A 

We then revert the above procedure using J \u\ 2 = 1 and 

£ 14 n^nwiic 1 - 



AC{1,...,N} keA 
N 



k4A 



N 



E E Xn{xj)Ylxn(xk)]j0--Xu(xk))=J2xu(x j ) 



j=l AC{1,...,N} 

which produces (|28p . 



□ 



Lemma 9 (Local uncertainty principle / Neumann Lieb-Thirring inequal- 
ity). Let u E S?oa &e an N-anyon wavefunction on M. 2 , and Q a square with 
area \Q\. Then 



[ 2N \Dju\ 2 XQ(xj)dx > C'AIqP)- 1 [ 



p(x) 



'Sol 
\Q\ 



dx. 



(29) 



Proof. We use \Dju\ > Vj|w| and then apply the Neumann Lieb-Thirring 
inequality given in Theorem 1 141 in the appendix to the bosonic wavefunction 
Id. □ 



Remark. Note that the limit case Q 

|2 



K 2 , N = 1, is just a special case of 
the Sobolev-type inequality HVdJij > C p ||u|| 2 4//( - p 2 ^||u||p P ^ p 2 ^ with p = 4. 

Lemma 10. For any domain ft C R 2 raf/i finite area \fl\ we have 



i 

P 2 



101 



2 



>(l-4e) / ^ + (2 _ £ -i ) (i|^)! ) (30) 



/or arbitrary e > 0. 
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Figure 3: Example of a splitting of Qo and a corresponding tree T of sub- 
squares. For the B-square at level 3 in the tree, the set A{Q) of associated 
A-squares consists of 8 elements, while for the two B-squares at level 2, A(Q) 
coincide and has 4 elements. 



Proof. First note that 
'SaP 



i 

P 2 



loi 



, 4 








>/ 




Jn 





7nP 

loi 



!-> 4 

2 



where the first integral on the r.h.s. is equal to 



i 

f)2 



JnP 

loi 



3 
P 2 



JnP 

loi 



+ 7 



101 



101 



Together with Holder's inequality applied to the negative terms 

-i i 

'SaP 



p 2 < 



2 , and 



P |0 



101 



this produces the inequality ([30 



□ 



Theorem 11 (Kinetic energy inequality for anyons). Let u G f 0a &e an 



N -any on wavefunction on R , u>i£/i N > 2. Then 

N 

V / |L>,«| 2 dx > Ca-C^ 
" Jm 2JV 



p(cc) 2 da;, 



(31) 



/or some positive constant Ck- 

Proof. Given u, and hence p, smooth and supported on some square Qo Q 
M 2 , we split up the domain Q := Qo into ever smaller subsquares according 
to the following algorithm (see Figure [3j): 
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• A given square Q is split into four subsquares Qj g {i 2 3 4} s -t- \Q'j\ = 
IQI/4- 

• Whenever Jq, p < 2, we will not split that square Q'j further, and we 
mark it A. 

• If all four squares Q'j are marked A, then we back up to the bigger 
square Q and mark it B. (Thus, we then stop splitting Q and the 
subsquares Qj are discarded.) 

• For each of the unmarked squares Q'j (i.e. those for which Jq, p > 2) 
we iterate the splitting algorithm with Q := Q'-. 

In the end (note that the procedure will eventually stop since p is integrable), 
we will have a system of subsquares Q organized in a tree T, such that: 

• The whole domain Qq is covered by squares of type A or B. 

• Let us denote, for a given B-square Qb, 



A(Q 



B) 



All A-squares Qa G T that can be found by going 
back in the tree from Qb (possibly all the way to 
k Qo) and then one step forward. 



then every A-square Qa G T satisfies Qa G A(Qb) for some B-square 
6 T. (Note that at least one B-square can be found among the 
leaves in the highest level of every branch of the tree.) 

• < Jq a p < 2 for every A-square Qa- 

• 2 < Jq b p < 8 for every B-square Q b ■ 

Finally, we also divide the A-squares into a subclass Ai for those subsquares 
on which p is approximately constant, and A2 for those subsquares with a 
non-constant density: 

Ai := jail A-squares Q G T s.t. f Q p 2 < c ( -^f- 

A 2 := jail A-squares QeT s.t. f Q p 2 > c-^f- 

for some fixed constant c > 1, to be chosen below. We will then split the 
kinetic energy integral T in the l.h.s. of (|31[) into a sum over all the marked 
subsquares Q G A\ U Ai U B forming the leaves of the tree T, and consider 
the three different types of squares separately. 
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Let us first consider the kinetic energy Tq b on any B-square Q = Qb- 
We split it up into two parts using k G (0, 1) and apply both Lemma[9] (local 
uncertainty principle) and Lemma [8] (local Pauli principle) to conclude 



T Q B > K 



C' 



> 



IqPJQ 
kC 2 



Sol 
\Q\ 



+ (!-«) 



c Q C l,N 



(1-46) | Q P 2 +^(2- e - 1 ) 



\Q\ 

+ (!-«) 



7c Q C ( 2 



a,N 



04 



cm; 

IQI 



where we also used that x — 1 > fS;^ 2 for 2 < a; < 8, as well as Lemma [TU1 



Choosing k and e appropriately, we conclude 



Tq b >CI n ^ Cl ^ 



B \Qb\ 



(32) 



with ci, C2 > 0. 

For the A2-squares Q = Qa, we use that the energy given by the local 
uncertainty principle is large due to a sufficiently non-constant p: 



'-Qa 



> 



IqP 



1 

pi 



' $QP 

\Q\ 



> 



a 



il - 46 - (6- 1 - 2K 1 ) I p\ 

JQ 



for < e < 1/4, where we again used Lemmas [91 and [TUl Taking e := 1/8 
and c := 24 we find 



Tq a > 



a 



p 



Q, 



(33) 



Lastly, we show that the remaining squares, of type Ai, have a negligible 
contribution to the energy compared to that already obtained from the Pauli 
energy on the B-squares. Namely, as observed above, every subsquare of type 
Ai is contained in A%(Qb) '■= <A(Qb) H A\ for some B-square Qb- Now, 
note that for each B-square £ T, say at a level k £ N in the tree, we 
have 

4 



(Iqb p) 2 



> 



\Qb\ -4- fc |Q |' 

while the total integral of p 2 over all Ai-squares associated with Qb is at 
most 



QeAi(Q B ) 



j=l Q€Ai(Q B ) 

at level j 



' \Q\ 



< 



E 3c 



4-J'|Qo| 



< 96 



4 fc+i 
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Hence, by ([32]), 

T Qb > Cl N Ll p 2 + | £ / P A . (34) 

Summing everything up, it follows from (|33|) and (|34|) that the total 
kinetic energy is 



d.r 



T = E /j^ n i 2 ( E xqA^)+ E xqbK)) 

i=i yR \Qa£^ Q B eB J 

> E + E > c xc\ n I 

Qa£A 2 Q B eB J ®° 

for a positive constant C^: := minjci, C2/96, C2/8}. □ 

Corollary (Lieb-Thirring inequality for anyons). Let u S 6e an iV- 

anyon wavefunction and V a real-valued potential on M 2 . Then 

N 

V/ {\D ]U \ 2 + V{ Xj )\u\ 2 ) dx > -C LT C- 2 N / |y_(x)| 2 da;, (35) 

~[ JR 2N ' JR 2 

for some positive constant Clt- 
Proof. The l.h.s. is bounded below by 

CkC\ n I p 2 dx- [ \V4pdx>--C K l C- 2 N [ V 2 dx, 

JR 2 JR 2 4 JR 2 

where we used f \V-\p < (fV 2 ) 2 (f p 2 ) 2 and minimized w.r.t. f p 2 . □ 
Remark. Theorem 1111 immediately implies the rough bound 

T > C K Cl N J Q P 2 dx> ( ^^ L (J n P- ldx ) = C « C In m (36) 

for the ground state energy of a non-interacting gas of anyons supported on 
a domain (ICR 2 . 



Appendix A: Improvements of the local energy 

Here we give some comments on how the explicit bound for the energy in 
Lemma [7] could be improved. One alternative approach for taking the full 
many-particle domain Q N into account is to extend the Hardy inequality (|18|) 
using a variant of Theorem 2.2 in [1], which here again has been modified 
to account for the underlying symmetry. 
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Figure 4: The set £1r, which can be mapped conformally to the unit disk 
by the map F. 



Theorem 12 (Many-anyon Hardy on a disk). Let VL := B\(Q) be a disk of 
radius A in M 2 and let u G @Q a . Then 



/ Y\ Dj u\ 2 dx>^£ Yf(?i±ZL, Xi ~ x i )\ u \ 2 dx, (37) 



where 
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f(R,r):= " (l + zJ^Cl-*) 1 7«_(i_ 2; ) 1+ T ii (i + 2 ) 1 
ThC 1 ~ R ) 

z is the complexification of the renormalized coordinate r with imaginary 
axis given by the unit vector R/R, and jr := 1 — - arcsin(i?). 

Proof. We proceed as in the proof of Theorem HI but this time parameterize 
not just Slo n, but the whole set by the relative coordinates R, r. For R 
at a distance XR, < R < 1, from the center of the disk, the total range 
of allowed r forms an eye-shaped set £Ir given by the intersection of two 
disks (see Figured]). After rescaling by (1 — i? 2 )~2A _1 and choosing axes 
so that the corners of the eye are at the points z = —1 resp. z = 1 in the 
complex plane, we can map the resulting set conformally onto the unit disk 
by the sequence of maps 

F : z i y — — — =: w, «)4 w~ =: C, C ^ t =: F(z), 

1 - z C+l w 

where jtt = jrtt is the angle at the corner of the eye. The resulting map is 
hence 

F(z) = y - ^ =: £ + i V , (38) 

(1 + z)i +(l-z)-y 

which is antisymmetric w.r.t. the antipodal map z \— > —z. Now, under the 
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coordinate transformation F, the integral w.r.t. r in (|2ip becomes 
/ \(D 3 - D k )u{r)\ 2 dr = [ KD^v^,^ 2 d^dr, 

>C 2 f ^ d£dn- ° l ' N f lF ' {z)l2 h\ 2 dr 
-^Js^e + V 2 ^- {l-R 2 )X 2 L \F{z)\ 2H dr ' 

where we again applied an annular decomposition of B\(0) and Lemma [21 
and used the antipodal symmetry of v(£,rj) := u(F -1 (£ + ir])). Finally, we 
have by the definition of / that \F'(z)\ 2 /\F(z)\ 2 = (1 - R 2 )f(R,r). □ 

Using the above Hardy potential / in (|26p . together with more precise 
estimates of the lowest eigenvalue of the corresponding operator H, and 
the restriction that the eigenfunction should be antipodal-symmetric, would 
almost certainly produce a significantly better bound on the constant 
and hence on the energy of a gas of anyons. 



Appendix B: A local Lieb-Thirring inequality with 
Neumann boundary conditions 



In this appendix we derive certain bosonic and fermionic kinetic energy 
inequalities on domains with Neumann boundary conditions. These follow 
straightforwardly from a recent method due to Rumin |18j . 

Theorem 13. Let Q be a cube in M. d with volume \Q\. Then 



N 



3=1 



p(x) 



'Sol 
\Q\ 



2(d+2) 
d 



dx, (39) 



where <f>j G H l (Q) and p(x) := Ylj=i \ 4>j( x )\ 2 - 

Moreover, if {4>j} are orthonormal in L 2 {Q), then 



N 



£l|v<^ 2(Q) > c d 

3=1 



p(*)*-\Q\- 



2(d+2) 
d 



dx. 



(40) 



Proof. For any e > and (f> £ L 2 (Q), let 

= ^+ + ^" e <+ := P { _ A ^> e} 0. 
and note that for <f> £ H l (Q) (interpreted in terms of quadratic forms) 



de 



P, 



{-A% >e} 

V V ' 

= f m l{\>e}dP(X) 



de < 



--f R \dP(X) 
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Denote the eigenvalues and orthonormal eigenfunctions of — Aq by {^k}kLo 
resp. {uk}^L . Then for each x G Q, 



\^-{x) 



} <t) (x) 



( E U k( x ) u k,' 
\X k <e 



E (u k ,4>)u k {x) 

\ k <e 
\A fe <e 



< 'l«l + ^<.l«l 



', (41) 



and hence \(j) e '~(x)\ 2 < + C^e^ ) ||0|| 2 by the well-known asymptotics 

for the Neumann eigenvalues on Q. 

Now, by the triangle inequality in C N we have for arbitrary {(f)j}f =1 C 



L\Q) 

N 



X 



> 



N \ 2 



a; 



, x e Q, 



+ 



and hence we find that for p(x) := J2j I0i( a; )| 2 an d {</>j}f=i ^ H l {Q) 

N !■ poo 

- (( Q 1 +C d ei) ! p{y) dy ) rh dx 
JQ 

i 



> 




> 



C d [ p(y)dy f f 
JQ JQ JO 

d 2 C d j Q P 



(d + 2)(d + A) J Q 



p{ x ) 



C d \Q\ 



1 \ 2 d 

— e4 



de da; 



1 \ 5 



C d |Q| 



2(d+2) 
d 



da; 



/?(a)2 - 



2(d+2) 
d 



das, (42) 



with C' d := d 2 C d d /(d + 2)(d + A). 
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In the case that {4>j} are orthonormal, Bessel's inequality applies in (|4ip : 



E 



u k {x)u k ,(t>j 



\X k <e 



2 




2 


< 


y~] uk(x)u k 






A fe <e 


A fe <e 



which then replaces every occurance of Jq p in (|42|) by 1. 



□ 



Remark (Generalization). If we know how Eo(e, a) := 5Z^A fc <e l^fcC 32 )! 2 be- 
haves for a given domain f2 C M d , then we can evaluate 



N 



n jo 



p(as) 



i(e,x) / p 



de da; (43) 



for either the Neumann or Dirichlet Laplacian on Q. Again, p in the r.h.s. 
is replaced by 1 in the case that {4>j} are orthonormal. 

Theorem 14 (Many-particle version). Let Q be a cube in M. d with volume 
\Q\, and let u G H l {W dN ) be an IS -particle wavef unction. Then 

i 2(d+2) 

( Sqp\ 



N 

~l Jm. dN Jq 



p(x) 



\Q\ 



dx, 



(44) 



where p(x) := YE =1 / Rd (jv-i) . . . , sc^-i, a;, Xj+i, . . . , a^)] 2 da; fc . 



Proof. We define for each x' = (a?i, . . . ,Xj-i,Xj+x, . . . ,asjv) £ 
collection of functions 

x i->- cj)j[x,x') := u[x\, . . . , Xj-i, x, Xj+x, . . . ,xn) 

in L 2 (Q), and proceed as in the proof of Theorem 1131 writing 



td(N-l) 



\V x (f)j(x, x')\ 2 dx 



for each x' G M. d( - N ^ , and similarly to ([H 

TV 

jd(iV-l) 3 

N 



b e .' + (x,x')\ 2 dx de, 



E 



faj,x')| 2 ds' 



<E / (\Q\~ X + C d ei) [ \4> 3 {y,x')\ 2 dydx' 
~[ Jm d ( N - 1 ) v J Jq 

= (|Qr X +C d ei) y 
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as well as using the triangle inequality on L 2 (R. d ^ N 1 - > ; C N ) 
' f 



l 

ix,x')\ 2 dx'\ > 



N 



3=1 



N 

|2 j / 



/ J2\^-(x,x' 



)| 2 cfe' 



for x £ Q. Hence, 

N 



> / / \V x (pj(x, x')\ 2 dx dx' 

> I T pix^-UQ^ + C.ei) ! p) 
JoJo V Jo J 



2 



de dx, 



IQ 

and dm then follows as in (H2D. □ 
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